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SDL  DATA  SERVICES 


INTRODUCTION 

This  report  lists  the  data  services,  computer  programs  and 
other  data  available  from  the  Seismic  Data  Laboratory  (SDL) .  The 
following  items  are  described  in  this  report: 

I.  Digital  computer  programs  (abstracts). 

II.  Digitized  seismic  data. 

III.  U.S.  Coast  and  Geodetic  Survey  epicentrrs 

on  punched  cards. 

IV.  Earthquake  bulletin  data  (LRSM  and  VELA 
observatories) , 

V.  Shot  Report  data  on  punched  cards  from 
82  U.S.  nuclear  explosions,  i 

Space  is  available  for  persons  who  wish  to  study  data  at  the  SDL. 
Those  who  wish  to  visit  the  SDL,  or  request  data,  should  direct  their 
inquiries  to: 

Headquarters ,  USAF/AFTAC 
VELA  Seismological  Center 
Washington,  D.C.  20333 

Att:  Project  Officer,  Seismic  Data  Lab 


!•  DIGITAL  CUMHulfcH  PROGRAM  ABSTRACTS 

ALL  PHOiJHmMS  i,j  IfilS  l  I S 1  HAVE  BEEN  JN  FORTRAN,  COPIES  OF  THE 

SOURCE  UtO*,  RKOjRAM  l  ISriNUS*  A ND  PRUGRaM  WRHe-UPS  ARE  AVAILABLE  ON 
REQUEST. 

A*  C'JHPfet.ATIUN  mud  HpgRESSIUN  ANALYSIS 


1,  AOtCfiV  AUTUCUvARIANCt  SUBROUTINE 

FnRlRAN.  COMPUTES  I  HE  AU  TUCJOV  *i«  I ANCE  OF  A  SERIES  OF  N  TERMS 
hi  IN  |.  LAGS, 

2*  AUfOCO*  AUTO-CORRELATION  ANALYSIS 

FORIRaN.  GIVEN  A  SERIES  OF  VALUES  X  C±  >  #  X (2 > # . . . X ( N ) ,  THIS 
PHUgRaH  ^ I Ll  COMPUTE  The  PKOuUcT-mOmenT  CORRELATION  COEFFICIENT 
(AUTo-cOHHelaT I  ON  >  BETWEEN  SUCCESSIVE  TERMS#  WHERE  THg  LAG*K> 
GOES  FROM  U#l*.*.*Lk 

3.  cRscnv  autu-  and  cross-variance  subroutine 

FOR  I  HAH •  GIVEN  TWO  SERIES  X  AND  Y  EACH  OF  N  POINTS*  COMPUTES 
AUfOCOVAKI ANCE  OF  X  ANU  Y  ANO  CRUSS-VaRIANCE  OF  X  TO  Y  AND  Y 
TO  X. 

4.  FRENORh  FHEUUENUY  AND  TEST  OF  NORMALITY 

FOR  I  RAN •  GIVEN  DIGITIZED  SEISMIC  noise  data*  this  routine 
INVEST l GA I tS  The  DISTRIBUTION  of  SEISMIC  NOISE  AMPLITUDES. 


LOGNORil 


MEAN  ARU  SAMPLE  VARIANCE  OF  POWER  SPECTRA 
ESTIMATES 


FOR  I  MAN .  To  COMPUTE  fHfe  SAMPLE  MEAN#  VARIANCE#  AND  CONFIDENCE 

limits  of  powep  spectral  estimates,  under  assumptions  that  loo 

OF  THE  POWER  IB  NORMALLY  DISTRIBUTED. 


manycor 


CORRELATION  ANALYSIS 


FORJHAN.  GIVEN  M  SETS  OF  DATA#  THIS  PROGRAM  CALCULATES  THE 
MEAN,  CORRECTED  SUM  OF  SUUAR6S#  VaRIarcE#  AND  STANDARD  DFVIaTION 
FOR  FACH  SfcT,  ALSU  CUMPUTfcS  IHg  REGRESSION  COEFFICIENTS# 
(LINFAR)  CUOREiaTION  COEFFICIENTS,  AND  THE  STANDARD  ERROR  OF 
FS I  1 Ma I E  FUR  ALL  POSSIBLE  COMBINATIONS  OF  SETS  OF  DaTa. 


7.  PULSCOK 


SI  A' istical  analysis 


FORIHaN.  To  CORHELaTe  A  GIVeN  PULSE  OF  A  SfISMOQRaM  WITH 
CORRFSPONUlNG  TIME  I  NT  fcH V  aLS  THROUGHOUT  THE  SEISMOGRAM  USING  THE 
DENOYEH  ME  I  HOD. 


a,  tfcucnR 


RECURSIVE  CORRELATION 


FORTRAN.  To  COMPUTE  I  HE  MINIMUM  VARIANCE  ESTIMATE  OF  THE 
JOINTLY  COHrfLaTFD  SIGNAL  FROM  MUl T I -CHANNEL  RECORDS. 

b  •  TIMF  SERIES  ANALYSIS 
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BUCK  Ya 


CNCJJS  Sf'feCI  ha,  A 'JO  POWER  bPfcCTRA 


fORlRA'**  This  PHU^RaM  IS  a  HfcVlSfcu  AND  SUPPLEMENTED  VERSION 
OF  T  Hfc  PfuiUHAM  Tl'KE  Y  F  ROM  UcSu.  C  APAH  I L I  T I  fcS  INCLUDE  FILTERING# 
DEC IMA I INu#  REMOVING  THeNq  of  DEOREt#  PLOTTING  OF  SeMILQG  aND/OR 
I  I  'EaH  Ptu's  Of  POWER  ANU/IJR  AUTOCORRELATION  SPEcT«a»  ANn  SPECIAL 
smooth i ng  feature, 

CUHERNCY  TIMfc  SERIES  ANALYSIS 

FORTRAm.  Tn  cOMPUTfc  ALL  THE  ORDINARY  COHERNC I ES,  a‘TTO* 

spectra,  and  phase  relations  of  a  set  of  input  data  channels. 

COlYTO*Y  SPECTRAL  ANALYSIS 

Fortran,  This  program  computes  auto-spectra»  cross-spectra, 
AUTo-CURRcLaTIonS,  CROSS-cURRELaTIONS  and  COHERENCIES  using  the 
COUL Y»T  UKt  Y  METHOD*  FOURIER  IRaNSFURmS  OF  THE  DATA  ARF  FIRST 
COMPUTED  and  Then  MANIPULATED  To  GlVfc  THE  DESIRED  RESULTS.  THIS 
PROGRAM  is  InTfNDFD  AS  an  ALTERNATIVE  to  PROGRAM  blacky. 

COOL  HYPtR-RAPlD  SPECIALIZED  COOLEY-TUKeY  FOURIER 

TRANSFORM 

FORTRAN.  To  cOMPUTfc  THE  FOURIER  SERIES  EXPANSION  OF  A  RFAL* 

or  complex-  vai ufd  data  series#  or  the  data  series  from  the 

COMPi EX-VALUED  FOURIER  SERIES  EXPANSION, 

COULER  HYPER-RAPID  SPECIALIZED  COOLSY-TUK£V  FOURIER 

TWANSF  ORM 

FORTRAN.  To  COMPUTE  The  FOURIfeR  SERIES  EXPANSION  OF  A  REAL 
VALUED  TIMfc  SERIES. 

COULTWO  FOURIER  TRANSFORM  OF  TWO  DATA  SERIES 

SIMULTANEOUSLY 

FORTRAN.  Tn  COMPUTE  THE  FOURIER  SERIES  EXPANSION#  USING  COOL 
(O.V.)  OF  TWO  nATA  SERIES  SIMULTANEOUSLY. 

fastfti  fast  direct  and  invehse  Fourier  transform 

FOR  I  RAN •  COMPUIES  I  HE  DIRECT  AND  INVERSE  FOURIER  TRANSFORM 
OF  A  PERIODIC  FUNCTION. 

FKSPTRUN  TIMfc  S£HI£S  ANALYSIS 

FORlHAN.  Tn  COMPUTE  AMD  DISPLAY  THE  FHFOUENcY-WaVE  NIJMREP 
PO*EP  SPEC'RA  np  SEISMIC  NUIS£  ALONG  WITH  A  RESPONSE  FUNCTION 
FUN  T Hfc  CuKhfSpOLDTnG  nh^Y.  iHlb  PROGRAM  IS  INTENDED  Ac  A 
nfcRLAC-MfcN 1  FOP  CERTAIN  PaRTS  UF  PROGRAM  PEAKAY  USING  Tmf 
C(l  JLr  Y-  |  Lftfc  Y  MpTHOu  10  ESTIMATE  aUTU  aND  CROSS  POWpR  SPFrTRAL 
'  e  m;  i  r  t  i  u,*ct  i nNS . 
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9,  FOUlAUhl  FOU«lfcH-LAUU£HKi;  TRANSFORM 

FORIRAN.  EXPANDS  A  GIVEN  T I  ME  FUNCTION  IN  A  SERIES  OF 
UUUcRrtfc  ►  Oimc  T  f  OMS*  ANU  ►RUM  THE  I.AGUeRRF  EXPANSION  COMPUTES 
FOUHIER  AfiML  I  TiiDE#  PHASE*  And  POWER  spectra. 

10.  F0ULAUR2  UIHtCl  ► OUHIEH-LAGUERRE  TRANSFORM 

EORIRaN.  GIVEN  A  TIME  SERIES.  THIS  PROGRAM  COMPUTES  AND 
PLUTS  ON  A  SEMJ-I.OG  OH  LUg-LOG  BASIS,  THE  AMPLITUDE.  POWFR, 

ANU  PHASE  SPECTRA.  The  PRUGRaM  PRINTS  THESE  QUANTITIES  ALONG 
WITh  The  I H  LAGUERPb  FUNcIlUN  EXPANSION.  AND  COMPUTES  THE 
WEIGHTING  ►UNCTION  IN  POSHjVE  OH  NEGATIVE  TIME  MEASURED  AT  THE 
LAUUeRNE  Sample  POINTS  of  The  GIVEN  TIME  SERIES* 

11.  FOUHaN  FOUHIER  ANALYSIS  SUBROUTINE 

forthan*  given  a  time  series  of  n  points,  computes  sine 

TRANSFORM.  COSjNF  TRANSFORM,  MODULUS  a.  modulus  normfd.  log  a, 

phase  (fraction  of  a  circle),  and  maximum  value  of  the  modulus. 

12.  FUURTR  FOURIER  ANALYSIS  SUBROUTINE 

►ORIRAN.  COMPUTES  • HE  SlNg  OR  COSINE,  SMOOTHEO  OR  UNSM00TH6D 
TRANSFORM  Up  A  SFRIES  OF  M  TERMS. 

13.  FOUTrS*T  FOURIER  TRANSFORM 

FORTRAN.  GfVFN  A  TlhE  SERIES  OF  N  POINTS,  THE  PROGRAM 
COMPUTES  AMPLITUDE,  PHASE,  FREQUENCY,  AND  A  POWER  SPECTRUM  WITH 
THE  USE  OF  FOURIER  TRANSFORMS. 

14.  framis  cohpuiaIion  of  spectra  AND  PHASE  from  The 

LAPLACE  TRANSFORM 

foriran.  given  the  parameters  in  the  laplace  transform,  this 
PHUGRam  evaluates  The  AMPLITUDE  SPECTRA  and  PHaSE  OF  the 

THmNSF'JHMA  I  ION. 

15.  FTP  HYPEH-RAPIO  FOURIER  TRaNSfORM  PaC«AGE 

foriran.  this  is  a  Package  of  six  subroutines  for  calculating 
FU’'RIEH  THanSFRRHS  and  MANIPULATING  COMPLEX  data.  TMc i R  NAMES 
are-  FIPACIV66,  DFTRNU63,  TAPEH63,  COOL,  SCALE63*  AND  DRUM63. 
the  package  allows  detrending  and  TaPeRInq  on  direct  transforms, 
and  Tapering  d4ta  on  inverse  transform,  in  admtion  the  package 
will  CUHPUlg  am  AMPLITUDE-PHASE  or  POLAR  REPRESENTATION  of  the 
FJUHIEH  TkanSFoRM  on  DIRfeC^  TRANSFORM  AND  WILL  ACCEPT  THF  OATA 
IN  DIFFERENT  F'nRMS  FOR  INVERSE  TRANSFORM. 

16.  FT^DcOOL,  ► T3DCOOL  TOO  AND  three  DIMENSIONAL  FOURIFR  TRANSFORM 

PACKAGE 

FORIRAN.  The  SUBROUTINES  IN  THIS  PACKAGE  COMPUTE  TWO  AND 
THnef  DIMLNSIUfiAL  FOURIER  TRANSFORMS.  THEIR  NAMES  ARE"  FTgDCOOL, 
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FT^DfOUL#  I'UQL#  MaTKaoj#  a^U  SCALE.  aS  WITH  CO0L»  THE  DIMENSIONS 
ON  The  da  ( A  mu«;T  bfc  A  POWER  01'  TWO. 

J.7 »  GHIZEL  FOUH'tR  TRANSFORM  BY  GOERTZE'L  S  METhOn 

F 0 R I R a N •  To  cOHPUTte  The  FOURIER  SERIES  EXPANSION  OF  A  REAL* 

valued  even  functiun#  ur  rHt  heal-valued  even  function  of  a  real- 
VALUfD  FOURIER  SERIES  EXPANSION, 


lb.  MULTICOH  TIME  beHlES  ANALYSIS 


FOR  f  HAN •  This  Program  COMPUTES  MULTIPLE  coherence  functions 
FOR  SEISMIC  ARRAY  DATA  RAPIDLY  AND  EFFICIENTLY.  GIVEN  AN 
ORIGINAL  StrT  OF  N  SUBSET  DATA  CHANNELS#  THE  PROGRAM  WILL  COMPUTE 
THE  N-i  MULTIPLE  COHERNCE  FUNCTIONS-  X < I > ( N- f /N, . . . , N- I ♦. ) , 

where  i  *  x . . . m— i .  this  program  will  then  reorder  the  n  data 
CHANNELS  any  number  OF  times#  EACH  Time  COMPUTING  ANOTHER  N-i 
MULTIPLE  COHERENCE  FUNCTION.  The  print-out  INCLUDES  A 
DESCRIPTION  OF  -The  NOTATION  used.  OPTIONAL  PRINT-OUT  INCLUDES 
ALL  THE  AUIO  AMD  CROSS  SPECTRA.  In  ADDITION  a  PROVISION  EXISTS 

tu  pi.o>  the  mUi  tjple  cuherence  functions,  the  cooley-tukey 
method  of  spectral  estimation  is  used  to  obtain  high  speed# 

19,  partlcoh  time  series  analysis 


Fortran.  This  program  computes  partial  coherence  functions 
f oh  taped  Data.  the  program  also  computes  the  amplitude  and 
phase  of  the  associated  transfer  function,  the  output  includes 
PRINTOUTS  OF  These  FUNCTIONS  as  well  as  plots  of  the  COHFRFNCE 
FUNCTIONS.  THF  COOLtY-TUKEY  METHOD  OF  SPECTRAI  ESTIMATION  IS 
USED  TU  CBIaIN  HIGH  SPEED, 


20.  PREDICT  TIME  SERIES  ANALYSIS 

FORiRAN.  To  cOMpUTt  A  PREDICTION  OF  THE  STRAIGHT  SUm  OF 
MULTI-CHANNEL  ARPAYS#  1SP&N  UNITS  AHEAD  IN  TIME  WHERE  ISPAN  IS 
THE  PRcD I C  *  ION  SPAN.  I  HE  PREDICTED  SUM  AND  THE  ACTUAL  SUM  ARE 

subtracted  to  yield  a  prediction  error,  and  all  are  plotted  to 
the  same  scale  factor,  also#  to  compute  and  plot  power-spectra 
of  these  i races. 


21.  PRER 


TlMt  SERIES  ANALYSIS 


FORTRAN.  THIS  PROGRAM  DOES  PREDICTION  ERROR  FILTERING  OF 
SEISMOGRAMS, 


22.  PSu 


THEORETICAL  POWER  SPECTRAL  DENSITY 


FORIRAN.  To  COMPUTE  THE  THEORETICAL  POWER  SPECTRAL  DFNSITY 

for  gaussian  markov  model. 


23.  RLlMPWR 


REAL  Time  POWER  SPECTRA 


FORTRAN.  To  COMPUTE  Tre  REAL  TIME  POW£R  SPecTRa  FUNCTIONS  FOR 
ElfHFR  uVEMLAPPING  OH  NON, OVERLAPPING  TIME  INTfRVALS  FOR  a 

specified  duration.  * 
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24.  TUKfc*.  TUKtf  SPECTRUM,  CROSS  SPECTRA  AND  POWFP 

SPEC'IRa 

f  or • ha  n  »  This  not  stales  analysis  program  contains  rwatF 

B  A  >  1 C  SObH^UTlNE?.  I  He  FI«S1  TWO#  FILTER  ANn  REMOVAL  Of  TREND# 
PRfcPAHt  lot  DATA  FOR  the  SPECTRUM  ANALYSIS  SUBPROGRAM.  TIME Y 
SPtCTHi'M  COMPUTES  FOH  iHt  TWO  SIMULTANEOUS  TlMP  SERIES#  TWf  CROSS 
<CU“  ANU  UUAnH»  TURE->  SPEC  *  RA  #  AND  THE  TWO  POWpR  SPEf'T‘A.  PHASE 
AND  COHEWt ARE  CALCULATED* 

25#  VELQLAG  LAGOtHRE  EXPANSION,  FOURIER  ANAL  IS  I S 

Fortran,  computes  ihe  laCueRRe  expansion  of  surface  waves 

RECORIJEU  M*  STATIONS  aLONQ  A  PROFILE.  FROM  THE  LAGUpRRf 

expansion  ii  computes  f uur i gR  spectra#  and  from  the  fouhieR 

spmctha  II  commutes  PHASE  VeIOCITT.  ATTENUATION#  and  0  AS 
FUNCTIONS  Op  F&EOUtNCT  bfcTWEEN  PAIRS  OF  STa>IOnS. 

26.  VF nSPThm  TIME  SERIES  ANALYSIS 

FORIRAN.  To  cUMPUTe  AND  DISPLaY  THE  FHEUUE*' CT-WA  Ve  NUMBER 
PO^EP  SPEC*Ra  ^F  SEISMIC  NOJSE  ALONG  WITH  A  RESPONSE  FUNCTjoN  FOR 
the  corresponding  vertical  ahhat. 

.  DIGITAL  FILLING 

1.  AWLGFLIR  ANALOG  FILltR  SUBROUTINE 

FORTRAN.  TwIS  SUbRUUljNE  PEHFURMS  LUR  PaSS  AND  HlQH  p4SS 
FILTFR1NG  Up  digital  data  In  THE  Same  way  that  analog  filtfps  DO 
OF  ANALOG  data, 

2.  CALC  PROGRAMMED  ARITHEMITIC,  HFAI 

FORTRAN.  GJVFN  IHfi  CENTER  FREQUENCY  AND  BANDWIDTH  OF  A 
BA^D-PaSS  MLTcR,  This  SUBROUTINE  COMPUTES  RECURSIVE  DIGITAL 
FILTER  COEFFICIENTS. 

3.  COtFFT  FILTER  COEFFICIENTS  FOR  BAND  PASS  FILTe« 

Fortran,  given  The  ctNi eh  of  band  pass#  half-rIdth  of  band# 

and  POLL-UFF  Iw  terms  Uf  NORMALIZED*  FREQUENCY  CFREOOEN rt  TIMES 

the  time  i ncpe*#ent i .  a^u  the  number  of  coefficients  in  half  of  a 
symmetric  p  ilTpr.  computes  and  returns  the  rand  pass  filter 

COtFFICIENfS. 

4.  DCUNVOLI  INVtRSE  CONVOLUTION  FILTER 

FORTRAN.  To  REMOVE  THfi  RESPONSE  IN  A  SEISMOGRAM  DUE  TQ  THE 
INHERENT  FILTER  CHARACTERISTICS  OF  the  SEISMOMETER,  FILTER, 
VELOCITY  TRANSDUCER#  AND  GALVANOMETER,  THE  ROUT  I NF  AMOVES  THE 
RESPONSE  UF  ONt  INSTRUMENT  at  A  TIME  USING  THfc  FILT£P 

chahac i eristic*  of  Thai  instrument,  sucessive  passes  ihpough 
THt  PO'illNt  n|ih  DlFFEPfcNT  FILTER  CHARACTERISTICS  WILL  PFMOVE 
ALL  THt  ImiESI »'Ep  RESPUNSES, 
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5.  F  A  'j  T  f  I  L  FASI  BROADBAND  MLTtH 

F  ( i  W  Ik  a  w  *  Tn  hHUAUBAiMli  PaSS  K  I  L  1  b  N  SE I  SPOOR A  MS . 

t>.  fCUN  GENERATE  FILJfch  COfcFFlfcMS 

FORTRAN.  This  ROUTINE  IS  UStL)  r U  GENERATE  COEFFICIENTS  FOH 
TH.;  filbuFu-le  LiGUEHHt  FILTER  <  P  U 1  N  T  F  L  T  >  SUBROUTINE. 

7.  F  1 L  L)  F  C  OIGITaL  FILTER  AND  DEC  I MATOP  SUBROUTINE 

FOR  I  RAN  #  FILTERS  AND  DBClMATbS  A  GIVEN  SERIES  OF  N  POINTS 
WUH  GIVEN  FILTER  CObFUCIbNfS  AND  RETURNS  THF  FILTERED  AND 
DEDlMAlfcD  SERIFS. 

»•  F1LNDIS  DIGITAL  FILTERING 


FOR  I  RAN#  Tn  FILTER  TlMt  SERIES  WHICH  ARE  READ  FROM  A  BINARY 
INPUT  rAPc  W I Tm  A  set  UF  FILTER  COEFFICIENTS  and  WRITE  THE 

filtfred  data  on  an  ouiput  binary  tape,  the  filter  coefficients 

MAY  PE  read  in  FROM  Ca«DS  UR  They  may  be  SYMMETRICAL  BAND-PASS 
COEFFICIENTS  DESIGNED  w 1 1 H  THE  USE  UF  SUBROUTINE  COfFFT  <UfS 
G1J2).  THE  PRnGPAM  has  THE  OPTIONS  OF  PRINTING  THE  FILTFR 
COEFF I C I EN I S , 


».  FILPI.OI  FREQUENCY  RESPONSE  OF  DIGITAL  FILTER 


FORTRAN.  l 
RESPONSE  UF  A 

coefficient. 


.  TO  COMPUTE.  PRINT#  AND  PLOT  THE  FHEOUENT Y 
SYMMETRIC  digital  filter.  AS  SPECIFIED  BY  ITS 


2.  ALTEHNA1 I VELY#  to  compute  the  filter 
COEFFICIENT  From  the  DESIRED  filter  CHARACTERISTIC,  given  the 
NUMBFR  OF  COEFFICIENTS#  AND  THEN  CARRY  OUT  <jj. 


10.  FILTER  filth  SUBROUTINE 


FORTRAN.  APPLIES  GIVEN  FILTER  CUEFFlClENTS  To  a  ScRlfS  OF  N 
POINTS  AND  RgTliRNS  THg  FILCHED  SE«lES. 

11.  landmi  a  recursive  filter 


FORTRAN.  r,IVtN  THE  COEFFICIENTS  OF  A  N U Mg R I (; a L  FILTER  AND 

consecutive  segments  OF  an  input  data  serifs,  THIS  ROUTIhc 
GENERATES  I  he  SEGMENTS  OF  The  CORRESPONulNG  F 1 1  T(  Rpp  SERIFS.  THF 
ENTIRE  SERIES  is  FILTERED  by  means  OF  SUCCESSIVE  calls  TO  THF 
SUBROUTINE.  3  '"f 


12.  L0PA7 


LOW  PASS  AND  HIGH  PASS  FILTERING  BY 

approximation  TO  analog  filtering 


DIGITAL 


FORTRAN.  TmIS 
FILTERING  up  Da Ta 
SUoSFT  TAPE  and  a 


program  Performs  low  pass  and  high  pas* 

ON  A  SOBSeT  Format  TAPE#  AND  PRODUCES  ANOTHER 
PLOT  I  APE  OF  THE  FILTERED  DATA. 
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13.  MAIN  INTEHPOLAT IONS  AND  APPROXIMATIONS# 

FILIbHING 

FORTRAN.  THIS  PROGRAM  PRODUCES  RECURSIVE  DIGITAL  FILTER 
COEFFICIENTS  FpR  LOW-PaSS  OR  BAND-PASS  FILTERS.  AS  SUCH  MOST 
OF  ITS  OPERATION  IS  I/O  AND  BOOKKEEPING#  WITH  THE  COMPUTATIONAL 
ASPECTS  OF  THE  PROBLEM  BEING  RELEGATED  TO  SUBROUTINE. 

IT  SHOULD  BP  NO  f  EU  (HAT  THfc  PROGRAM  DAS  WRITTEN  WITH  LOW-PASS 

AN'J  PANd-PhSS  INCORPORATED  subroutines  to  generate  digital 
FILTfRS,  however#  n  Can  be  used  To  model  a  GENERAL  FILTER 
provided  1*E  COEFFICIENTS  OF  TH£  CONVENTIONAL  PROTOTYPE  FILTFRS 

TRmNTFER  FUNCTION  ARE  GIVgN  AND  THE  POLES  Of  THE  PROTOTYPE  ARE 
OF  MULTIPLICITY  ONE. 

14.  MAXLIK  MAXIMUM  LIKELlHUOD  FILTER  PROGRAM 

fortran.  This  program  computes  and/or  applies  a  21  op  30 

POINT  HAX1HUM-I  IKLIHUOU  F JLTER  (REALIZABLE  UR  SYMMETRIC  >  TO 
SEISMIC  AH«aY  PATa.  PRJNTfcD  OUTPUT  IS  PRODUCED  AS  WELL  AS  A  PLOT 
Tape  DISPLAYING  the  input  and  OUTPUT  DATA  ANn  USEFUL  COMPUTED 
information.  in  AUDITION,  a  Save  Tape  IS  GENERATED  WITH  all 
PERTINENT  UATA  ABOUT  The  OPERATION  written  IM  AN  EASILY 
ACCESSIBLE  FOR*-..  A  PRUVISION  EXISTS  TO  RUN  THE  PROGRAM  FROM  THE 
SAVE  TARE  PkOCcSSING  NEW  DATA  WITH  BOTH  KINDS  OF  MAXIMUM- 

likelihood  filteps.  in  this  case#  eithep  filter  is  ffficifntly 
re-computed  from  information  written  on  the  save  tape  and  a  new 
SAVE  Tape  Is  GeNFRAIED  with  DETAILS  OF  the  NFW  FILTER  and  its 
OPcHATlON. 

The  PPIMTEO  OUTPUT  INCLUDES  the  actual  summation  CONSTRAINT 
SATISFIED  ay  Th6  FILTER  COEFFICIENTS.  THIS  IS  FOR  REAL ! 7ABLE 
FILTERS  AND  fOp  SYMMETRIC  FILTERS. 

15.  KAX-LIKE  MAXIMUM-LIKELIHOOD  FILTER  PACKAGE  66 

F  OR  *  Pan.  ThIS  IS  a  PACKAGE  OF  SUBROUTINES  F OR  COMPUTING  A 
MaX1MUm-LIKr;LIhOOU  F1UER  ►  OR  DIGITaL  DATa.  GIVEN  THE  CORRELATION 

mairix  of  a  noisf  Sample  and  a  summation  constraint  matrix,  only 

THfc  CALLING  SEQUENCE  T"E  FIRS!  aNd  MASTER  SUrROUTINe#  LEVIN66 

IS  DFSCPIBED,  CINCE  aLL  T he  OTHERS  IN  THE  PACKAGE  ARE  CALLED 

internally  fro*  it  or  from  each  other,  their  subroutine  names 

AHfc-  bLKF.Gc,  Maine#  MaIMPY#  move,  MAIDEN,  and  SIMEO7.  THE 
RESULTING  UlTpR,  when  APPLIED  ro  A  GROUP  OF  INPUT  CHANNFlS  OF 
DA  I A ,  WILL  PRODUCE  ONE  UR  MORE  OUTPUTS  SUCH  THAT  THE  TOTAL 
OUTPUT  POWER  IP  a  MINIMUM  SUBJECT  TO  THE  CONSTRAINT  SUPPLIED. 

this  filter  is  sometimes  referred  to  as  a  minimum  variance 
UNBIASED  FILTER.  LEVINSON  RECURSION  IS  USED  TO  COMPUTE  THE 
filtfr.  ihjs  procedure  Tares  advantage  of  The  symmetry  of  the 
matrices  in  order  to  increase  the  maximum  number  of  input 

CHANNELS  and  The  ALLOWABLE  FILTER  LENGTH  WHILE  DECREASING  THE 
RE'JUtRfcu  COMPUTING  TIME. 

16.  PHILTER  LINEAR  MULTICHANNEL  CONVOLUTION  FILTER  FOR 

TAPED  D*TA 

FORIHM.  Am  INPUT  (APE  OF  UNMOLT  I  PLEXfiD  DATa  IN  THE  SUBSET 
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17. 


t 


FOUMH  1  S  oOMV'LVfcU  W  |  I H  A  MUL I  1  CHANNEL  LINEAR  FILTER  AND  WRITTEN 
f )  Nj  i  u  A'  C'jlpnr  TAPt  IN  THE  SUbbfcT  FURMAT. 


!»• 


POINTFL  I 


DISCRETE  l.AGUERRE  FILTER 


FOM I  HA  'M  •  TwlS  SUbRUUIjNE  ACCbPlS  ONE  OAT*  INPUT  AT  A  TIME 
ANU  PASSES  Tm  I  <!  DATA  THROUGH  A  SET  OF  RECURSIVf  FILTER 
OPERA  HONS  WHlrH  CAN  Bt  STATED  EITHER  AS  A  TIME  OPERATION  OR  AS 
A  frequency  OPfRaTION. 


REOF  IL 


INTeHPOLaI IUN  AND  APPR  x I  mat  I  ON#  FILTfRINQ 


FOR  I  RAN  •  GTVFN  AN  1NPUF  TIME  SERIES#  A  CFNTER  FREQUENCY,  A 
BANDWIUTF'#  AND  THE  DATA  SAMPLING  RATE  IN  SECONDS  PER  POINT#  THIS 
SUBROUTINE  HpTiiRmS  The  FILTERED  time  SERIES#  The  FILTER 
COEFFICIENT#  t>  NO  A  COMPUTED  PARAMETER  OF  THE  FILTER  POWER 
FUNCTION.  the  program  will  also  filter  data  in  SECTIONS# 


IV.  PECFIL3  SINbLfc  PASS  PHASE  SHIFT  RECURSIVE  FILTER 

FORIRAN.  Tn  PERFORM  UnE-PaSS  PHASE  SHIFT  RECURSIVE  FILTERING 
OF  ScjShIC  IJATa.  RIVEN  I  Hfc  (J-VaLUE#  CENTER  FREQUENCY#  AND  SAMPLING 
RAIE. 

20#  R&SRPC  PLOl  RESPONSE  OF  RECURSIVE  RAND  FILTER 

FORIRAN.  COMPUTES  AND  PLOTS  THE  AMPLITUDE  OF  A  RECURSIVE 
PHASPLPSS  UIRITAI  FILTER. 


21.  TWA  SPAUAl  INIERPOLATION 

FORTRAN.  TP  COMPUTE  A  LeaST-MEAN-SOUARe-eRROR  FILTER  WHICH 
I N I  ERPUI.A  I  PS  OhE  CHANNEL  OF  AN  ARRAY  FROM  CERTAIN  OTHER  CHANNELS 
IN  THE  ARRAY,  APPLY  the  filter  to  CONSTRUCT  an  interpolated 
fstimaIe#  uonstriict  the  error  time  function,  and  display  the 

PO«*ER  SPEC  I  HUM  OF  THE  tRHOR. 


D.  data  RETRIEVAL  ANn  maNJPULaUOn 

1.  APCTTAl  REARRANGE  subset  data 

FORIRAN.  This  PROGRAM  obtains  all  possible  COMBINATIONS,  in 
SE1S  OF  Two,  OF  SUBSET  DATA  CHANNELS  AND  FORMS  AN  OUTPUT  TAPE 

containing  a  label  with  each  set  of  two  records# 

2.  CENTRE  cenieh  a  record  subroutine 

FORTRAN.  TwlS  SUBROUTINE  TENTERS  THE  RECORD  Ain  of  n 
POINTS  IMU  THf  PECORD  B<I>  OF  M  POINTS  AND  PLACES  IEROES  ON 
BOTH  SIDES. 

3.  CNTUR4  CONIOURING 

FORTHAN.  CnTURA  WRITES  PRINTER  PLOT  OF  THE  CONTOURS  OF  A 

two-dihensionai  array,  the  contours  are  formed  by  shading  the 

area  BtTWbfcN  LEVFL  2*  AND  LEVEL  21  ♦  j,  THIS  DISPLAY  IS 
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SUPfcRlHROSto  UPON  ElTHfcH  A  GRID  WORK  OR  A  BORDER  WITH  TIcK-maW* i, 
THc  rONTOUH  lEVELS  HAY  BE  SPECIFIED  INDIVIDUALLY#  OR  BY  THE 
INCREMENT  BETWEEN  LEVELS.  IN  ADDITION*  PROVISION  HAS  BEEN  Made 
TO  At  LOW  I  he  ORIENTAI  ION  OF  THE  ARHaY  TO  BE  ALTERED  BY  El  THE11 
E  *,Jh  A  Nu  i  NO  THE  ROWS  ANU  COLUMNS  UR  REVERSING  EITHER  The  POW‘>  OR 
THE  OULIJMNS, 

4,  COLLATE  MERGES  EPlCENlbR  LISTS 

For  Iran*  This  program  Iakes  two  binary  Tapes  or  a  racked 
B  HapY  TAPE  ANn  EPICENTER  CAROS  and  searches  FOR  EVENTS  WHICH 
ARfc  RERORIED  on  roth.  matching  is  determined  by  comparison  with 
(i>  A  PREDETERMINED  ORIGIN  time  DIFFERENCE  WHICH  INCLUDES  A 
DEPTH-TIME  CORRECTION  and  <2>  A  PREDETERMINED  DISTANCE  DH'L«ENCE 
which  includes  a  latitude-longitude  correction,  latitddf. 
luugitude#  depth,  magnitude,  and  origin  time  differences 
(ReSidUaLS*  Apr*  COMPUTED  ONLY  F  OR  MATCHING  PaRTS  OF  EVENTS  AND 
AHt  AVAILABLE  AS  OUTPUI  . 

5.  DEPThmAG  RETRIEVAL  AND  COUNT  OF  EPICENTERS  SUmmaRI7ED 

by  depth  and  magnitude 

Fortran#  This  hrugham  has  the  ability  to  retrievf  data  from 
ant  onc  of  four  differently  formatted  tapes,  the  ability  to  save 

ALL  CARD  images  RETRIEVED.  ABILITY  TO  SPECIFY  RETRIEVAL  CRITERIA 
RELATIVE  ru.  A.  HECORUING  stations 

B.  DATA  SOURCE  codes 

C.  DISIANCE  and  AZIMUTH 

D.  OTHER  RELATED  MISCELLANEOUS  FEATURES 

6.  DETHnD  REMUVAL  OF  MEAN  AND/OR  LINEAR  TREND 

FORTRAN#  THIS  SUBROUTINE  REMOVES  THE  MEAN  ONLY  OP  Ht  Th  'HE 
MEAN  AND  LINEAR  TREND  UF  A  SERIES  OF  N  POINTS,  AMO  RFTl'R'S  THf 
DETRENDED  SERlrS. 

7,  eplist  epioenter  listing 

FORTRAN.  LIST  the  INFORMATION  CONTAINED  IN  pee -punched 
USCGS  bP  I CENTEP  CA^OS  FQH  EARTHQUAKES  SATISFYING 
PREDETERMINED  WAI  UbS  FUH  ANY  DESIRED  COMBINATION  OF  THE 
FOLLfiWlNG  VARIABLES-  LATITUDE.  LONGITUDE,  DFPTH,  TT'ipF , 

AZIMUTH,  ANJ  DISTANCE. 

B.  MERGSEIS  MeHGINg  OF  i  SEISMOGRAMS  UMdT.-R  ONg  LABEL 

Fortran.  This  program  will  take  any  two  seismograms  from  one 
TAPE  and  MERGE  THEM  TOGETHER  UNDER  A  NEW  LABEL  WITH  50  CHANNEL 
IOcNTIL IFRS. 

9,  NORMAL  NORMALIZE  SUBROUTINE 

FORTRAN.  COMPUTES  ANU  H£TUHNS  a  NORMALIZED  SERIES  FROM  A 
GIV.W  SERIES  OF  N  TERMS, 
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10.  «A  JGF  maximum  MINIMUM  VALUE  SUBROUTINE 

FORTRAN.  FlNOS  AND  RETURNS  THfc  MAXIMUM  AND  MINIMUM  VALUES  OF 
A  SERIES  OF  N  TERMS, 

11.  SQUASH  SUALINQ  ROUTINE 

FORTRAN*  FINDS  AND  RETURNS  THE  MAXIMUM  AND  MINIMUM  VaL.UeS  OF 
A  SERIES  OF  N  TERMS  aNU  Then  SCALES  THE  SERIES  WITH  A  G I VgM 
MAXIMUM  ABSOLUTE  VALUE. 

12.  TAREaCI  EPICENTER  DATA  STANDARDIZATION 

Fortran.  Various  epicenter  cards  from  several  sources  are 
WRITTEN  onto  tapf  in  a  si aNDARU  format  IN  BCD  OR  PACKED  binary 
MODE.  Sfc  1 S  ANn  GEO  REGION  NUMBERS  MaY  BE  ADDED  INTO  THE  FORMAT. 
BINARY  OR  BCD  TAPES  MAY  BE  LISTED  AT  END  OF  AN  UPDATE  RUM  OR 
MAY  be  CONVFRTfD  BACK  IU  STANDARD  FORTRAN  BCD. 

13.  TAPMFRG2  TAPE  M£RGE 

Fortran*  To. merge  iwo  magnetic  tapes,  made  up  of  epicenter 

card  IMAGES  on  MAGNETIC  TAPE»  THE  RESULTING  MERGED  TAPg  IS  IN 
CHKONOLUGICaL  SEOUENCE  and  in  BCU  mode  WITH  END-OF-FlLPS  BETWEEN 
months  and  A  DOURLE-ENU-0F-FIL6  AT  THE  END  OF  THE  YEAR. 

14.  TRAVEL  I  CUMPUTER  P-ARRIVAL  TIMES 

fokiran.  computes  ihe  p-arrival  times  from  seismic  events# 
REPORTED  on  pdc  cards  uf  THE  u.s.  coast  and  GEODETIC  survey,  to 
specified  stations  and  lists  them  together  with  all  the 

INFORMATION  cOmTaINED  on  the  pde  cards. 


*  INTERPOLATION  and  APPROXima t ions 


1.  aUIOdEU 


con i hulled  solution  of  ordinary  differential 

6UUAT IONS 


For  Iran,  to  provide  a  convenient  means  uf  solving  a  set  of 
N  SIMULTANEOUS  FIRST  ORUfcR  DIFFERENTIAL  EQUATIONS,  The  RaSIC 
integrator  is  a  self  starting  variant  of  the  adams  method 
INCORPUF'A  HNG  AUTOMATIC  STEP  SIZE  CONTROL,  USFR  PROVIDED 
STATEMENTS  are  EXECUTED  AT  SPECIFIED  VALUES  OF  EITHER  THF 
INUEPENubNl  variable  oh  any  CONTINUOUS  FUNCTION  of  the  DFPFNDENT 
variables,  automatic  control  uf  truncation  or  discontinuity 
error  is  provided,  the  dependent  variables  are  internally 
carried  in  partial  double  precision  to  control  round  off  epror. 

2.  RES  BESSEL  FUNCTION 


fortran. 
J(N>  OH  I ( N ) 
N  IS  A  FIXED 


This  routine  evaluates  for  a  given  x  and  n,  either 
where  x  is  a  normalized  floating  point  mumbfr  and 

PnlMT  INTEGER. 
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3.  HFFGh 


SHhfcNlC*L  BESSEL  FUNCTION 


FoRIRaN*  This  R0UI1NE  CALCULA  TfcS  THE  spherical  bessel 

FUNCTION  0*-  THr  FIRST  K  I  N[)  AND  OF  T  HE  SECOND  KIND  TOGETHER  WITH 
THe  F I  T  AND  PEC0ND  DERIVATIVES  F  uH  a  GIVEN  APGUMENT, 

4.  CERES  ELEMENTARY  FUNCTIONS 

F'OHfRAN.  TuIS  SUdHUUlINL-  WAS  WHITTEN  TO  PROVIDE  THE  SINE* 

COS  I mE»  A  HU  EXPONENTIAL  FUNCTIONS  PRECISELY  ACCURATE  TO  H 
SJGNiFlCANl  DECIMAL  FIGURES  for  ARGUMENTS  no  GREATER  THAN  ROUGHLY 
4U.i4> , 

5.  CLABMU  EIGENVALUES  aND  EIGENVECTORS 

FORTRAN.  3l VEN  A  CUMPlbX  SUUARE  MaTHIX#  NOT  NECESSARILY 

hehmitian.  to  find  all  eigenvalues  and  one  eigenvector  for  each 

DISTINCT  EIGENVALUE.  » Hb  ROUTINE  USES  The  La  PUDDE  METHOD  WITH 

de term i nan i  evaluation  and  the  mueller  method. 

6.  CilPTPWW  POLYNOMIALS*  EVALUATION  OF  A  SPECIAL  FUNCTION 

FORlRAN.  This  IS  A  FUnCTIUN  SUBROUTINE  WHICH  EVALUATES  I  HE 
EXPHFSSION  for  THE  POWER  Al  A  GIVEN  FREQUENCY  for  A  -Band  pass 
FILTFR. 

7.  6XTREM  MULT I-plIENSIUNAL  EXTREMUM  SEEKER 

FORlRAN.  COMPUTES  I Hfc  MAXIMA,  MINIMA,  OR  SADDLE  POINTS  FOR  A 
SINGLE  FUNCTION  OF  UP  10  gU  PARAMETERS. 

8.  GLUUAU  GAUSSIAN  LEGENDRE  QUADRATURE 

FORTRAN,  This  ROUTINE  IS  DESIGNED  TU  OFFER  A  SUBSTANTIAL 
IMPRoVEMENi  0 VpR  THE  USUAL  USE  OF  SIMPSON  S  RULE  TO  APPROXIMATE 
DEFINITE  INTEGRALS,  THE  GAUSSIaN  LEGENURE  METHOD  HaS  THE  HIGHEST 
ALGEBRAIC  PRECISION  FOR  A  UlVfcN  NUMBER#  N,  OF  EVALUATIONS  OF  THE 

iniegrand. 

9.  GRSCH  VEC'UR  ORThO nORMALIZaT  JON 

foriran.  This  subRuut i ivg  orthonormalizes  The  input  set  of 

VECTORS  by  MEANS  OF  THE  GRAM-SCHMIJT  ORTHOGONAL  I ZaTION  PROCEDURE. 
WHILE  UOING  THIS.  IT  CMEOKS  EaUH  NEW  VECTOR  COMPUTED  TO  SFF  IF 
ITS  NORM  IS  lE^S  THAN  A  SePcIFIEU  QUANTITY.  IF  IT  FINDS  A  NULL 
-VECTOR,  AUU  HEmCF  A  LINEARLY  DEPENUfcNT  INPUT  SFT,  IT  RETURNS  AN 
ehhor  flag  anu  terminates  IhE  calculation, 

10.  iufum  interpolated  function  and  its  derivative  AMD 

INTEGRAL 

Fortran,  givfn  any  stT  of  abcissas  in  monotomcally 
INCREASING  UPDcR  with  correspond! ng  ordinates-  The  INTERPOLATED 
functions  are  OBTAINED  a!  I  HE  DESIRED  ABCISSAS-  OR  THE 
INlEiPULA IfcD  OPR  I Va  T 1 VfcS  QF  TRt  FUNCTION-  OP  THE  INTEGRALS  OF  THE 
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FUNCTIONS-  fJR  THF  SUM  UF  Thk  1  (m  fb&^AL  •  Tut  INTEGRALS  OF  THE 
FUNCTIONS  and  THE  SUM  UF  TheSc  INTtCiHAlS  AH'-  Al  WaYS  COMPUTED. 

the  curve  is  Parabolic  between  adjacent  pairs  of  points  amp  t m 

D  E  R  I  V  a  I  I  V  t  ANO  FUNCTION  IS  MATCHED  ON  BOTH  SIDES  OF  EACH  POINT, 

11.  LAPHIV  A  SUBROUTINE  FOR  FINDING  THF  LAPLACE  INVERSE 

OF  A  COMPLEX  FUNCTION 

FOR  I  HA'N  •  G I  VpN  A  FUNCTION  Of  The  COMPLEX  VaRIABLF  S,  aND  a* 
ARRAY  OF  I  he  ImDPPENDENT  VARIABLE  time,  to  GENERATE’  THE  TIME 

function  fob  each  of  the  time  values  indicated. 

12.  LSUPOL  LEASi’  SUUaReS  POLYNOMIAL  FITTING  SUBROUTINE 

FORTRAN.  GIVEN  X < p , X ( i ) , X ( 2  > » • • . » X ( M )  AND  F  <  I  > ,  F  <  2  ) »  . . . » 
F(H>,  WHERE  F  C I  >  IS  The  OgSp.RVED  DEPENDENT  VARIABLE  AND  X<I>  IS 
THE  OBSERVED  DcPFNDENT  VARIABLE  ANU  X ( I )  IS  THF  OBSERVED 
INJEPENDENI  VARIABLE#  THE  POLYNOMIAL  Y  »  B<1>  +  p<2>  *X  ♦  t.. 

+  BU1)*x**X  IF  FITTED  FOR  ALL  DEGREES  OF  K,  FROM  K«i  TO  k  « 
K(MAX)  W  i  I H  CFRTAIN  OPIlONSi 

13.  RFH  PROGRAMMED  ARITHMETIC,  REAL 

fortran,  Tpis  subruijiine  employs  the  regijla  falsi  method  in 
FIND  the  Value  of  THE  INDEPENDENT  Variable  for  a  PARTICULAR  VAI  UF 
OF  THE  DEPENDENT  VARIABLE  (OR  FUNCTION). 

1*.  RKAM  NUMERICAL  SOLUTION  OF  ORDINARY  DIFFERENTIAL 

equations 

FORTRAN.  To  COMPUTE  iRfc  NUMERICAL  SOLUTIONS  OF  A  SYSTEM  OF 
N»  (N.LE.lug),  SIMULTANEOUS  FIRST  -ORpER  ORDINARY  D I  FFE  RFNT  I  Al. 
EQUATIONS  WITH  GIVEN  INIIIaL  OONUI r ions  over  a  given  INTFRVAl 
using  partial  double-precision  arithmetic,  thf  method  of 
SGLUTIUn  may  bp  EITHER,  <!>  RUNGE-KUTTA,  (2)  adams-moulton. 
WIlHOUl  ERROR  CHECKING,  OR  (3)  AUAHS-MOULTON  WITH  error  CHFCK1MG. 
THE  LATTER  TWO  USE  RUNUfc-KUTT A  STaRTErS. 

15.  RKAMSUB  NUMERICAL  solution  OF  ORDINARY  DIFFERENTIAL 

EQUATIONS 

FORTRAN.  TO  COMPUTE  IhE  NUMERICAL  SOLUTION  OF  A  SYSTFM  OF  N, 
TN  * LE • 10 U ) ,  SlMUL TaNEUUS  FIRST -ORDER  ORDINARY  D I FFERE NT  1 Al S  w|Tm 

given  initial  conditions,  using  partial  doubi e-precision 
arithmetic. 


SPEEnUR 


VfeClORS  AND  SIMULTANEOUS  LI  'Ear  Equations 


FORIRAN.  to  DETERMINE  AN  N 


MAXIMUM  ABSOLUTE  VALUE  UF 
P  -  AX.  here  p  is  an  u  x 
Rank  n  wpeRe  m.le.n  ♦  lt 
MAX  or  cFEBYSHpV  SOLUTION 


vecTUR  X  Which  MiNjMi7gS  The 


all  the  components  in  the  vprTnp 
1  VECTOR  and  a  is  an  m  x  n  matrix  of 
an  x  which  o°ts  this  is  callfd  a  mjn- 
10  the  OVER-DETfcRMINFD  SYSTpM  AX  ■  R. 


PAGE  12 


17. 


THAP7 


fKAHL/OlDAl.  H*  ’L  h:  INTEGRATION 


FORIRAN.  EVALUATES  THfc  IflltfjHAL  OF  P  <  X  )  U  X  FOR  A  SUCCF  SS  I  OR 
OF  M  X )  S  S*  THe  THAPEZOIOAl.  PULE. 

18,  ZOtPpftiT  ZUEPPH I  I Z  S  fcOUAflONb 

FORIRFN.  To  CALCU|.A)b  COEFFICIENTS  UF  RfcF  LFC  T  l  ON  AND 
RFF'RaCUuw  fop  GIVEN  BOUNDARY  CONDITIONS  USING  ZOFPr>H  I  T  ’/  r» 
EQUAT I u  r  j  S  • 

F,  GENERAL  SEISMIC  PROGRAMS 

1.  ARESPONS  OFF"btANElJ  ARRAY  RESPONSE  PROGRAM 

FORTRAN*  TmIR  PRUGHaH  COMPUTES  aND  PLOTS  Tug  UB  LOSS  Fl'A 
SIGNALS  C.UM  I NG  IN  FROM  a^gLeR  DIFFERENT  FROM  Tm£  |jjpAMSTf.|  p 
AZIMUTh.  ELABORATE  PRINTOUTS  ANU  PLOTS  are  GENERATED  IN  OPIIFR 
TO  make  The  output  SELF-EXPLANATORT . 

2.  CAL  IR  MAGNIFICATION  CALCULATION  USING  UlGlllZHj 

SINE  WAVE  CALIBRATIONS 

FOR  I  RAN •  COMPUTES  IHE  AVERAGE  PtAK-1 O-PeaK  AMPLITUDE.  THF 
MAGNIFICAT »ON.  AND  NORMAL  I Z I NG  FACTORS  UF  A  DIGITIZED  SlHf  Wa  VF 
CAL IrtRA  T 1  UN ,  THF-  Y  FACTORS  for  EACH  TRACE  MA>  BE  INPUT  OR  THE  Y 

may  re  calculated,  given  specified  INPUI  PARAMETERS. 

3.  C0N2A  CONVOLUTION 

FORIRAN.  To  ado  two  SPIKE  SEISMOGRAMS  1U(,fTHF-R.  CONVOLVE 
THEM  WITH  a  THIRD  SEISMOGRAM  ANr  CONVOLVE  Thf  RESULTS  wjTh  A 

wavelet  TO  produce  a  simulated  seismogram, 

4.  DEGOSTE  GHOST  reflections 

FORTRAN.  To  REMOVE  GHOST  REFLECTIONS  WHEN  a  SURFACE  AMD  DFFP 
WELL  TRACE  AT  THE  SAME  GAJN  IS  AVAILABLE. 

5.  dE'.Sfis  Vertical  aRRaY  processor  package  ii 

FORIRAN.  TO  06GHUSI  VERTICAL  ARRAY  SF J SMOG* A  MS  GIVEN  an 
ECHO-TIME  AND  PEFLECHUN  COEFFICIENT  FOR  EACH  SEISMOnPam.  JN 
ADuITIUn.  a  RUv  TRACE  UF  THE  UEEPWELL  TRACES,  AND  A  RIJH  traCC  Of 
THE  aLIGNEU-DEGHOSTEU-UEEP-WELL  TRACES  ARE  COMPUTED.  AL^O  THE 
COMPnNtNT  WHICH  IS  JUINTLY  cURReLaIED  On  aLL  Al  IGNED-DF:GwoSTfD- 
DEbP-WtLL  IRACFS  IS  COMPUTED.  FURTHERMORE  THIS  PROGRAM  Pl.OTS  ALL 
DEEPWELL  AND  DFGHOSTEli  JEgPwELL  TRACES,  THE  PAIR  OF  SlJM  TPaCFS, 
ANU  THE  CORRELATION  TRACE, 

6.  DISCPE I E  EVALUATION  UF  kaUTe  EXPANSION  COEFFICIENTS 

FORIRAN.  EXPANDS  EITHER  A  SglSrilC  SIGNAL  OR  A  SPECIFIED 
FUNCTION  IN  TERMS  OF  A  SET  OF  URThONORMAL I  ZED  EXPONENTIAL 
FUNCTIONS.  PARAMETERS  GIVEN  BY  THE  USER  DEFINE  THfc  OrThonoRMAL 
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K Oi«C T 1  ON-S •  THE  SPECTRUM  THfc  EXPANSION  IS  AVAILABLE  AS  A N 
op i low. 

7.  D|UPlA»  DISPLAY  AND  ANALYSIS  OF  AVOMALIES 

FORTRAN*  To  TAKE  THfc  OUTPUT  FROM  I  HE  ANOMAl.Y  PROGRAM  AND 
COMPUTE  THfc  AVERAGE*  SUN^ard  DEVIATION*  AND  NUMBER  OF 
OCCURRENCES  OF  A  SfcLtC'EU  ANOMALY  FOR  SELECTED  EUBARRAYS  OR 
SENSORS,  u  i  STANCES*  ANU  AZIMUTHS* 

t».  DISTA  UlSUNCt*  A/IMITH,  AND  TRAVFL  TIME 

FORTRAN.  COMPUTES  I  HE  OISIANCE  IN  DEGREES  AND  KlLOMETFRS, 

THE  AZIMUTH,  Back  AZIMUTH,  the  TRAVEL  TIME  OR  ARRIVAL  T I MF ,  AND 
THr.  ROUTIONS  PAPAMET£HS  FROM  A  SPECIFIED  SEISMIC  EVFNT  TO  A 
GIVEN  ST  A  I  ION. 

9.  DSST  DUMMY  St  1 SMUGR A  m  SUBSET  TAPE 

FORTRAN.  CREATES  A  DUMMY  SUySfc T  TAPE*  REMOVES  MEAN  UF  FACH 
CHANNEL  ANU  NORMALIZES  ENERGY.  IHfcN  SMOUThS  EACH  with  a  filter. 

10.  OESPlKt.  0ESP1KE  SEISMOGRAMS 

FORlRAN*  TO  PF.MOVE  SPIK[-S  FROM  StlSMOGHAMS  BY  SIMPLY 
I NSEP I  1 nG  A  COP  I  ME  FUNCTION  IN  A  SPECIFIED  INTERVAL* 

11.  FL1PT  ELL  1 P 1  I C I T  Y 

FORTRAN.  To  rUMPUTt  I  HE  N„iIU  UF  THE  RADIAL  AND  VERTICAL 
COMPONENTS  OF  RAYLEIGH  MOTION  AS  A  FUNCTION  OF  PERIOD. 

12.  GRUVpL  GRUUp  VELOCITY  FROM  PHASfc  VfLOcITY  DIVERSION 

Fortran,  this  program  Takes  the  phase  velocity  versus  pehiod 
OH  FREUUENCY  POINTS  ANU  USES  Thfi  POLYNOMIAL  AND  ITS  DE«P'aTiVfS 
TO  OpTaIN  The  GROUP  VfLUUITY  a s  A  FUNCTION  OF  PERIOD  or 
FREQUENCY. 

13.  HARKRIUR  RAYLEIGH  WAVE  DISPERSION 

FORTRAN*  To  COMPUTE  HAYt.ElGH  WAVE  GROUP  VELOCITY,  PHASE 
VELOCITY,  Mixn  inversion  coefficients  AS  A  function  OF  frequency. 

1 A .  HEFAI.UHP  MEAOUHEU  NUlSE  ISOTRUPIC  PROCESSOR 

FORlRAN.  TmIS  PRUgRaM  COMPUTES  ANDZUR  APPLIES  A  MULTICHANNEL 
iSulPUPlC  PROCESSOR  TO  SE I  Si Ic  ARRAY  DATA,  AN  ACTUaI  NOT  SF  MODEL 
IS  USED  COciPUTru  FROM  The  Spectra  of  a  SPecIFIFD  data  Sample. 
EITHER  A  Pup.iT  UR  A  CISC  SIGNAL  MODEL  CAN  BE  COMPUTED.  THE 
prugpah  then  solves  thc  multichannel  wieneR-hopp  equation  in  the 
FREQUENCY  domain  TO  ge1  The  OPTIMUM  FILTER  which  REJECTS  the 
noise  and  Passes  the  signal,  the  filter  is  written  on  a  save 
Tape  fur  future  IISe*  an  OPTION  EXISTS  To  FILTfR  a  GIVFN  PIECE 
OF  data  ANU  plot  the  FILTERED  TRACE  AND  THE  DIRECT  SUM.  PERSON 
PLUTS  ARE  ALSU  GENERATED  AND  PLOTTED  WHICH,  IN  CONJECTlON  WITH 
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iHt  OUTPUI  PLOTS#  ALLOW  iHfc  USER  T(J  CALCULATE  SIGNAL  TO  WOISF 
RA I InS. 

1*>.  ISUFIL  THEURtTlOAL  ISOTROPIC  PROCESSOR 

FORTRAN.  This  PROGRAM  COMPUTES  AND/OR  APPLIES  A  MULTICHANNEL 
ISOTROPIC  PROCpSSOH  To  SEISMIC  ARRAY  data.  AN  ANNULAR  ring  noisf 
MODEL  AND  CITHER  A  PU1NT  OR  A  UlSc  SIGNAL  MODEI  CAN  R£  SPECIFIED. 
THc  PROGRAM  THcN  SOLVES  I  Rfc  MULTICHANNEL  WIEN£-P-HOPF  FOUaTION  IM 
THt  FflhOUeNcY  nOHAIN  TO  G£T  THE  OPTIMUM  FILTER  WHICH  REJFcTS  The 

noise  amo  passes  the  signal.  The  filter  is  written  on  a  save 
tape  for  future  use.  an  option  exists  To  filter  a  given  piece  of 
DAlA  AND  PLOT  THF  FILTERED  trace  and  the  dINEC1  SUM.  PERSON 
PLUTs  mRE  also  GFNEHaTeo  AnO  plotted  WHICH#  IN  CONJUNCTION  with 
The  oUfPUl  PLOTS#  ALLOW  T Ht  USER  To  CALCULATE  signal  TO  NOISF 
RATIOS. 

16.  LAGTlME  Signal  and  NOISE  ADDITION 

Fortran*  This  phugRam  combines  a  signal  and  a  noise  onto  the 
same  data  ghanuei  at  varying  signaL-to-noise  ratios,  a  tape  "f 

THE  new  COM0jN4riON  UF  DATA  IS  WRITTEN  AND  A  PI  OT  TAPE  IS 

fuhmfd. 

17.  LOCATE  DETERMINATION  OF  LATITUDE,  LONGITUDE#  DFPTH, 

AND  ORIGIN  MME  OF  A  SEISMIC  SOURCE 

Fortran,  calculates#  s*  an  interactive  least-squares 

PRUCFSS,  I  me  LaUIUq;...  LONGITUDE#  DEPTH,  aND  ORIGIN  TIME  OF  A 
SEISMIC  SUURCE#  AND  DETERMINES  THE  CONFIDENCE  INTERVALS  and 
JOINT  CONFlOENrE  REGION  F OF  THE  SOURCE  COORDINATES.  COMPRESS¬ 
ION!  WAVE  ARRIVAL  TIMES  F^OM  F  I  Vfc  OR  MORE  STATIONS  ARE  MSFI)  • 

1«.  HAAAPAY  ARRAY  SIMULATION 

FoRlRAN*  PPODOCLS  A  SEI  of  OUTPUT  TRACES  CORRESPONDING  To 
an*  SEl  OF  ELEMENT  POSITIONS  DESIRED.  A  NOISE  BACKGROUND  IS 
SIMULATED  BY  SfLFCTING  TEN  TRACES  FROM  A  SET  OF  TWENTY  SaMPLF 

traces#  anu  assigning  io  each  trace  an  arbitrary  velocity, 

aZIMIIIH  anu  PElATIVfc  AMPLITUDE.  A  SIGNAL  IS  SELECTED  FROM  a 
SEl#  US  VtLOClTY  ANU  AZIMUTH  ARE  CHOSEN  aND  A  8 1 GNaI- TO-NO I SE 
RA  I  In  IS  assigned. 


jy.  MAaVOSb 

fohIhan. 

ESI Ima ' ES. 


mean  and  Sample  variance  of  spectral  estimates 
to  compute  the  f-ean  and  sample  Variance  of  spectral 


?u.  MOu  PV-7  SUHFACE  have  dispersion  and  amplitude 

FORTRAN.  This  PRUUHaM  lomputes  for  all  modes  OF  l  OVt  AND 
PaYLfIGH  WAVES  ON  an  ELASTIC  HALFSPACE  OF  PLANFrPARALLEL, 
HOnOGEf'EOUS,  I*oTRUPlc  LAYERS,  THE  FOLLOWING-  PHA^E  VfLOcITY# 
GROUP  VfcLUujTY,  AND  SURFACE.  ORBIT  aS  FUNCTIONS  OF  PERIOD  OR 
FREQUENCY-  AMPLITUDE,  the  PRODUCT  OF  VERTICAL  WAVE  NlJMRFP  TIMES 

l a tEr  thickness,  anu  potentials  as  functions  or  depth-  average 
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KHfcTH*  POTENTIAL.  AND  I  0  1  aL  HNfcH'iY  ULNKITIeS,  AND  AVERAGE 
H  U  '<  l"  7  0 1  1  /iL  gNEPUV  FLUX#  CORRESPOND  I NG  TO  EACH  LAYfcR-  AND  A 
SU‘'hA  I  ;  Ol'i  Of  THE  fc'VfcHG  T  UUANUUfcS  t  ROM  THE  F  RUE  SURFACE# 

RKOV | S i UN  IS  MaDF  FOR  PRiNlfcH  PLOTS  OF  PHASfc  AMD  GROUP  VELOCITY. 
PROVISION  iS  Al  sn  MADE  PUR  PARTIAL  DECOUPLING  mheRE  AMPLITUDE  is 
1.  a HGF-  JO  CERTAIN  CHAMRtLS  at  DEPTH#  RELATIVE  TO  THE  SURFACE 
AMPLITUDE. 


j>j  .  NE  1  W D H E  2 


EVENT  PROCESSOR 


KiHlRAN*  GiVFN  SETS  Uf  SEISMIC  STATIONS  AND  EVENTS#  THIS 

prjgRa''  processes  information,  relative  to  various  specified 
event  I'aGNITHDfS  and  a  threshold  probability,  this  PROGRAM 
rooPuTis  i he  network  probability  and  station  capability  of 
ntlErTJON  UR  I nENT IF  I C  A  T I  ON  OF  SEISMIC  EVENTS.  IT  ALSO  COMPUTES 

the^ threshold  magnitude  for  all  stations  and  UP  TO  !7  EVENTS  and 
subset*.  OF  events#  in  audition#  it  computes  The  number  of  events 
occurring  above  A  given  Magnitude  and  how  many  of  those  fvfnts 
will  tu  identified  oh  go  undetected  fop  specified  subsets  of 

EVENTS  aNU  STATIONS. 

22.  NOSPIKI  SPIKE  REMOVER 

FOR  I  RAN •  Tn  REMOVE  SPIKES  AND  CLIPS  FROM  SEISMOGRAMS. 


23.  PSuaI  PH  PSEUitU  DEPTH.  ALPHA#  AND  BETA 

FORTRAN#  TO  COMPUTE  THfc  PSfcUUO-LAYER  THICKNESSES.  PSFUDO- 
COiIRhESS  I  UNaI  AND  "SHEAP  VELOCITIES  fop  A  GIVEN  EARTH  STRUCTURE# 
ANY  IaYEH  mITHIN  THE  STRUCTURE  MAY  BE  DIVIDED  INTO  N  SUB-LAYERS. 
THfc  CORRESPONDING  ALPHAS#  B£TaS#  THICKNESSES,  DEPTHS#  AND 
DENSITIES  ARE  DEFINED  FOR  EaCH  SUB"LaYER  USING  A  LINEap  GRADIENT 
ACROSS  THE  ORIOINIAL  LAYER. 


2^.  respond  array  azimuth  response 

FOR i Ran*  To  cOMPUTfc  I  HE  RESPONSE  OF  AN  A«PaY  TO  WAVE  NUMBERS 
ALONG  ONE  AZIMUTH. 


25.  RESPONSE  AHRAY  RESPONSE 

FDRfRAN.  Tn  nALQULATfc  AND  PLUT  IN  K-SPACE  THE  RESPONSE  OF  AN 

array  of  seismometers  with  respect  to  the  response  of  tHe  same 
array  to  WAVES  with  INFINITE  PHASE  VELOCITY  <EX.  NUMBER  of  db 
DOWN  FROM  AERO). 


26.  RMSOML  '•  COMPUTATION  of  RMS 

FORTRAN.  To  READ  UP  To  9  SEISMOGRAMS  FROM  A  MAGNETIC  TAPE# 
IN  SUBSET  Tape  FORMAT,  DETREND#  FILTER  (OPTIONAL)#  COMPUTE  RMS 

on  a  particular  desired  region#  and  the  mean  of  the  rms. 


27.  rotate  SEISMUgHAM  POTATION  PROGRAM 

fortran.  This  phugnam  takes  a  subset  of  seismograms  with  the 
z#  r#  i#  and  Possibly  time  channels  and  makes  one  uf  three 
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POaSlMLt  <uT A  (  l  OoS  ON  I  Ht  PRtufc  TERM  1  NEU  SEISMOGRAMS  ON  TME  SUBSET 
TAPE.  IHfc  TrHfE  POSSIBLE  ROTATIONS  ARE- 

1.  U>  AX1«?  POTAllUN  GIVING  an  OUTPUT  Z,  Hit  AND  IT Z 
l.  U-'  PRIME)  AXIS  H  ( j  f  A I  1 U  N  GIVING  AS  OUTPUT  ZT#  RZT,  AND  TZ 
3.  </-H  PRIME*  AXIS  MUTATION  GIVING  AS  OUTPUT  ZW#  RZ#  AMD  T ZR 

OUTPUT  N  ILL  IUCLUlJt  a  NiiW  SUBSfcf  TAPE  WITH  TMF  Z,  R,  AND  T 
CHANNELS  <"r  POUR#  IP  *  He  flMfc  CHANNEL  IS  GIVEN)#  AND  A  PRINTOUT 
OF  SlIlU'feT  INFORM/ T  ION  AND  PLOT  SCALE  FACTORS. 

26.  SUtG^fcl  VtHllCA*-  AHMaT  PROCESSOR  PACKAGE  II#  AD  JUSTFD 

POHIRAn#  To  REMUVb  GHOST  REFLECTION  ON  VERTICAL  ARRaY 
SfelSROi.RAMS  aLTGHED  FOR  FIRST  P-MUTION  GIVEN  ONLY  Thf  FCPO-TIMES 
aNu  Ppt  L  t  C  T ION  COEFFICIENTS,  IN  AUUITIUN  SUM  TRaCFS  OF  pOTH  THE 
AL1GMEU  Dt-tPHEl  L  AND  ALluNItU-URGHUSlbD-DEEPHFLi  TRACES  A^F 
CU.IPilTEc#  mLSh  T Ht  COMPONENT  WHICH  IS  JOINTLY  CORRELATED  ON  ALL 
Dfc*#HOS  <  tD  StjSMOf.HAMS  IS  c"MPUTeu.  FUHTHbHMORF  PLOTS  ARF 
OB  I A I  I't*  D  UP  aLi  pEbPRbLL  AuO  UEGHOiTtD  UFEPwFLL  TRACES#  TMF 
SU^FnC*  THacF#  The  TWO  aHOVr;  SUM  TRACES,  and  The  CORRELATION 
trace  . 

2V.  SEE  SPEC  1 P 1 C  DRIVING  PROGRAM  FOR  F H AM  1 S 

POMRAH*  Tm  1  f  PROGRmI!  ‘  aS  wRIIIEN  FOR  THE  SPECIFIC  PURPOSE 
UF  DP1vING  FRA^IS.  PRaMIS  EVALUATES  THE  AMPLITUDE  AND  PHASE  OF  A 
tr *N*sp i  h  puMrTrop  given  ihe  parameters  Of  that  function,  thf 

INPUT  1(J  Stfc  the  UPEHa'ING  PARAMETERS  FUR  FRAM1S  AND  THE 

RARAMFthRS  OF  THE  LARLaCE  TRANSFORM.  THE  OUTPUT  OF  SEE  TS  A 
IISTIW  UF  amP|  1TUDE  CRaW  «H  N0RMaL17eD  TO  A  REOUfcSTFD  F^EOUENCY) 

A Nu  PH#  SI  AS  A  FUNCTION  UF  TRtOUENCY.  CASES  Ma  V  Bb  RaTCMEP. 

3t).  SEGMENT h  SEGMENT  MULTIPLEX  TWO 

FnRlRAN.  LONG  SblSMUGH/iMS  ARE  BROKEN  UP  INTO  MULTIPLEXED 
SbGMF  e  T  S  UF  MEM  FEU  LF.NGTh. 

31.  SlMUFEP  VtMlCAL  ARRAY  SEISMOGRAM  SIMULATION 

FORTRAN*  To  SIMULA  1 1  T»'b  UttPWELL  TRACE  OR  TRACES  FROM 
MEaSUhF  MENlS  Or  a  SURFACE  nr  i,e  aR-SURFACE  SEISMOGRAM,  GIVEN  Thf 
FChU-T  i  MtS  jjpT  w£F  N  THE  Depf  well  AND  THE  SURFACE  AND  HI  V h i  Thf 
REFLECT  1UN  COEFFICIENT.  T  *'F  SUM  TRACE  OF  THE  DEEP  WELL 
SEISPOMRAMS  and  TP.t  SUM  I  RACE  OF  1  HE  SIMULATED  DFtPWELL 
Sfc  1  SmUIRAMS  aRp  rUMPUTEl).  PLOTS  ArtE  OBTAINED  «'F  The  DEEouFLI 

setgponrams.  Simulated  dee^well  seismograms,  the. surface  tracf# 

ANU  the  two  ABOVE  SUM  IRARPS. 

3?.  Sl'*Gfj  SIHULAUD  GHOST  REFLECTION 

FORTRAN.  To  SIMULATE  NORMAL  DEbPWELL  TRACE*'.. 

33.  SI-TVaRM  VERUCAl  AHHaY  SIMULATION 

Fortran,  in  simulate  at  n  compunent  vertic/l  apray  according 
TU  The  FOLLokUG  CRITERIA- 

1#  SURFACE  TRACE  CUNT  a  )  NS  a  SIGNAL  COMPOSED  OF  a  PSEUDO 
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p-HUl  St  ANU  WAVEl-fcT  FORMED  BY  SURFACE:  RE VEPBEKA T 1 ONS  AT  SOURCE* 

2.  EACH  UPfcPWbLL  IHACt  MAS  SAME  SIGNAL  AS  SURFACE  WITH  THE 
addition  ue  surface  ghost  at  the  vertical  array* 

3.  EACH  OFEPWELL  SIGNAL  IS  SHIFTED  BY  THE  DF  SI  RED  TlMP  DELAY 
IN  UHDfcR  IU  SIMULATE  WHaT  WOUU  BE  RECORDED  IN  NATURE. 

4.  NO  UTHER  PULSES  ARE  CONSIDERED. 

5.  NOISE  IS  90  PER  OEM*  COHERENT  BETWEEN  SUCCESSIVE  ARRAY 

elements. 

34,  SPECAVE  SPECTRAL  MATRIX  AVERAGING  PROGRAM 

FORTRAN.  This  program  computes  and  averages  spectral  matrices 

FO'<  MUL  riCHAfyNpL  SEISMIC  AhRaY  DATA.  THE  OUTPUT  OF  ThIS  PROGRAM 

is  a  tape  in  Tme  same  format  as  the  program  -hpFaLUmp-  save  Tape, 
using  that  TaPf  As  input  to  the  -hefalump-  program  allows  the 
user  TO  cOMPUTp  MULTICHANNEL  FILTERS  FROM  an  ENSEMBLE  of  noise 
spectrm.  this  program  is  also  designed  to  add  to  or  update  an 
ENSEMBLE  SPECTpaI  matrix  ON  Tape  from  a  PREVIOUS  run. 

Zb*  SPfcCMAl  THEURET1CAL  purer  SPECTRAL  DENSITY  MATRIX 

F'  OR  IRAN.  To  COMPUTE  HE  THEORETICAL  POWER  spectral  DENSITY 
maikix  for  the  gaussun  mahkov  model,  all  spectra  are  normalized 
to  The  apsouitp  maximum  value  of  the  spectral  set.  in  addition 

EACH  SPECTRAL  COMPONENT  P<HS.f  ‘  OF  THE  MATRIX  IS  SCALED  RY  THE 
FACTOR  (L0GF(1.  +  IMHS.F  ><1>  *100.  >  >/lPTR^»FH1)  SO  THAT  AT  HIGHER 
FREQUENCIES  THF  SPECTRAL  VALUES  ARE  OBSERVABLE  ON  THE  MATRIX 
PLOT.  THE  POWfR  SPECTRAL  DENSITY  MATRIX  CONSISTS  OF  THE  REAL 
PART  BELOW,  and  THE  AUIO-SHHCIRA  along  The  diagonal. 

36.  SPtCMTHX  THEORETICAL  POWER  SPECTRAL  DENSITY,  COHERENCY 

AND  PHASE  MATRICES 

FOR  I  PAN  *  Tti  cOMPUlt  THEfiRpTlHAL  POwER  SPECTRAL  UfcNSlTjfS, 
COHERENCIES,  and  PHASES  F  pH  The  GAUSSIAN  MARKOV  CHAIN.  TWO 
MATRICES  ARE  GtNFRATED'  HAV1MG  aLL  POWER  SPECTRA  NORMaLI7ED  TO 
THE  MAXIMUM  v At UF  OF  THE  SPFCTHaL  SET,  COMERFNrY  TO  A  PLAUSIBLE 
VmLUf  OF  l »  ANn  PHASm  TO  A*  APPROPRIATE  VALUp  OF  PI.  IN  ADDITION 

THE  power  spectra  pirs.fi  is  scaled  according  to  the  factor 

CLOQF  «!•  ♦1P(PS.F>1*1U0.)/1P‘RS»F > 1 >  FOR  OBSERVABILITY  AT  HIGHER 
frequencies. 

37.  SPECSUHg  FIL'EH  AND  SUBSET  PROGRAM 

FORIRAN.  Tn  pETPEND.  FlLThR  (FOUR  POLE  BUTTERWORTh > ,  AND 
SL'BSfT  DATA  fRom  EITHER  A  LIBRARY  OH  SUBSET  TaPE  AND  PLOT  THIS 
Da  I A  •  QPUONS  ape  INCLUDED  TO  DEMAGNIFV  THE  DATA.  GENERATE  PLOTS 

of  The  results,  and  process  mure  than  one  seismogram. 

38.  VERPROC  VERIICAL  ARRAY  PROCESSING 

FORTRAN.  Tn  PROCESS  VgHTICAL  ARRAY  DATA  IN  THE  FOLLOWING 
MANNER- 

1-  dethcnd.  taper,  and  corhect  for  instrument  gain, 

2- cumpuIe  r.s  noise*  signal  peak  to  pfak  fxoUrsjon,  and  sin 

RATIO  FUR  ALL  W a  ELEMENT  S . 
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.5-  l-OHM  THE  DEGHOSTINU  OPERATOR  Gtjf  OF  W  FOR  EACH  VA  TRACE 

wuh/withuut  brennfr  weighing. 

4-  COMPUTE  THE  FOURIfcR  SERIES  EXPANSION  GT(J>  OF  LaMHA  OF  EACH 
G<J>  OF  W, 

•>-  CONVOLVE  GT(J>  of  LaMDA  WITH  its  corresponding  DFEPWFLI 
trace  to  PRODUCE  a  DfcGHOSTbD  CONVOLUTION  RECORD* 

xmng  spectral  matrix  estimates 

FORTRAN.  This  IS  a  PACKAGE  OF  THRfcE  FORTRAN  63  SUgROUT I NfS 
FOR  COMPOTIng  AN  ESTIMATE  OF  THE  SPECTRAL  MATRIX  FOR  N  CHANNELS 
OF  TAPtD  DATA.  THE  NAMES  OF  THE  THREE  ROUTINES  IN  THE  PACKAGE 
AHfc-  *MNU»  cOOLFH#  ANU  SMOOTH.  IN  ADDITION  TO  THRSF#  Three  MORE 
SUdHoUTjNES  aRf  aSSUMEU  To  BE  ON  THu  SYSTfcM  TAPE*  They  A«F-  COOL. 
DISC63*  ANU  FRASF. 

YFaCT  T  factor 

FORTRAN.  TO  COMPUTE  Y  CALIBRATION  FACTORS  F  0R  DIGIT! 7£D 
SElSMOdHAMb, 


II.  Digitized  Seismic  Data 


The  SDL  library  of  digital  data  includes  short  and  long  period 
seismograms  from  over  100  U.S.  nuclear  explosions  and  1000  earth¬ 
quakes  recorded  at  various  LRSM  sites  and  the  VELA  observatories. 

A  complete  list  of  digitized  data  is  available  upon  request,  but 
users  are  encouraged  to  visit  the  SDL  to  ascertain  that  the  data 
is  applicable  to  their  project. 

Digitized  data  can  be  punched  on  cards  or  written  on  tape  in 
a  7  track  IBM  compatible  format. 

III.  USC&GS  Epicenters 


Information  contained  in  the  USC&GS  PDE  cards  has  been  punched 
on  cards  and  written  on  magnetic  tape.  A  copy  of  the  tape  containing 
epicenter  information  from  1960  to  date  is  available  upon  request. 
FORTRAN  programs  to  retrieve  data  from  the  tapes  include: 

DEPTHMAG  -  retrieves  epicenter  data  which  meet  various 
criteria  such  as  location,  depth,  magnitude 
distance,  azimuth,  geographic  region,  and 
recording  stations. 

ARRIVAL  -  arrival  times  of  23  seismic  phases  at  any 
given  station  for  any  set  of  epicenters. 

TRAVEL  T  -  computes  arrival  time,  distance,  and  azimuth 
for  a  given  station. 


IV.  earthquake  Bulletin  Data 


Earthquake  bulletin  data  from  the  LRSM  teams  (February  1962 
to  date)  and  the  VELA  observatories  (February  1963  to  date)  has 
been  stored  on  magnetic  tape.  The  data  includes  the  PDE  cards  plus 
all  phase  information  for  those  stations  which  recorded  the  earth¬ 
quake.  Recorded  phase  arrivals  not  associated  with  an  epicenter 
are  also  on  .the  tape. 


V.  Shot  Report  Data 

Shot  report  data  from  82  U.S.  nuclear  explosions  have  been 
punched  on  cards.  These  reports  contain  phase  arrivals  recorded  ai 
the  LRSM  stations  and  VELA  observatories.  The  information  recorded 
includes : 

-  Shot  Name 

-  Location 

-  Magnitude 

-  Origin  Time 

-  Recording  Station 

-  Phase 

-  Arrival  Time 

-  Amplitude 

-  Period 


